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Abstract 

We calculate the resistivity contribution of tilted grain boundaries with varying parameters in sub-lOnm diameter metal 
nanowires. The results have been obtained with the Boltzmann transport equation and Fermi’s golden rule, retrieving correct state- 
dependent relaxation times. The standard approximation schemes for the relaxation times are shown to fail when grain boundary 
tilt is considered. Grain boundaries tilted under the same angle or randomly tilted induce a resistivity decrease. 

1. Introduction 

A widely used model to identify the resistivity contributions due to different electron scattering mechanisms in metal 
interconnects is the Mayadas-Shatzkes (MS) model Q. In addition to a bulk term, this model includes partially diffusive 
scattering at the wire boundary surface and scattering at grain boundaries (GBs), the two dominant scattering mechanisms in 
small diameter wires 0-0. The MS model can be used to fit measured resistivity with different diameters: a certain diffusive 
boundary scattering probability, GB density and barrier height such that the relative contribution of surface scattering and GB 
scattering can be obtained. Regardless the simplicity of the model, it has many shortcomings that become worse when the 
wire diameter is reduced. Eirst, the Euchs-parameter is used and represents the probability for diffusive scattering at the wire 
surface & 0- This parameter is not directly based on physical properties of the boundaries, e.g. does not take into account 
the roughness of the wire surface and neglects the quantum-mechanical aspects of confined electron states. Secondly, the 
GBs are represented as potential Dirac delta barrier planes perpendicular to the transport direction, which limits the scattering 
possibilities to Umklapp scattering. The resistivity is expected to be different when other orientations are considered. Einally, 
the electron relaxation times (RTs) are obtained approximately, neglecting the coupling of different electron states and their 
RT. 

In 0 was shown recently that the resistivity contribution of GBs in thin metal nanowires (below lOnm diameter) is 
much larger than the surface roughness contribution, while they are found to be almost equally responsible for the resistivity 
contribution in larger diameter nanowires. The correct state-dependent RTs are obtained and surface roughness scattering 
probabilities are retrieved from Ando’s model 0, describing surface scattering quantum-mechanically with use of the roughness 
properties of the boundary surface 0. 

We extend this work here by considering tilted GB planes and obtain the correct electron RTs and corresponding resistivity. 
We consider two models for tilted GB planes: tilted GB planes having the same fixed angle and randomly tilted GB planes 
without correlations between them. 


IT Model 


The conduction electrons are modeled in the same way as in 0, as particles with effective mass approximation (effective 
mass denoted by to*) in a box with square cross-section (diameter D). The box is an infinite potential well with periodic 
boundary conditions in the transport direction, for which the wire dimension is much larger than the transverse dimensions. 
The Boltzmann transport equation for a stationary state and the resulting expression for the current are given by: 
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with Ez the electric field (along the transport direction), kz the wave vector along the transport direction, the sub-band 
indices for the transverse dimensions and fny.,ny (kz) the distribution function for a specific sub-band. 






The first order deviation from the Fermi-dirac equilibrium distribution using with zero temperature assumed, is given 
by: 
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with = ±Y fcp — (ttUx/D)'^ — {irny/D)'^ the positive and negative wave vectors at Fermi level for each sub-band and 

'^he corresponding RTs. When the scattering probabilities are determined by Fermi’s golden rule, the RTs are coupled 
and have to satisfy the following equations: 
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with V the potential representing the deviations of the nanowire from the ideally conducting box. The sub-band state with 
reverse wave vector —has the same RT as the state with kz due to reflection symmetry of the equations. This implies that 
when only perpendicular GBs are considered, which limits scattering to the Umklapp process, that the equations for different 
RTs decouple. Assuming that the RTs are equal when solving Q, an assumption that one often uses in simulations (e.g. in 
pO)), gives the correct result (see Fig. (a)). However, when tilted GB planes are considered, other scattering events are 
allowed and the equal RT assumption does not hold (see Fig. 0 (b)). 

The RT solutions of these equations determine the conductivity a and resistivity p =1!a through: 
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It is clear that any way of approximating the sub-band dependent RT solutions has a direct impact on the resulting resistivity 
and its scaling. 

In what follows, we only consider GBs to be present in 1^ = Vgb, while being modeled as Dirac delta barrier planes 
distributed along the wire: 
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with Ugb, Lgb determining the scattering strength, representing the energy barrier and width of the GB respectively, while the 
functions Za{x,y) determine the positions of the GBs. 

In |[T|, there are N perpendicularly oriented GB planes distributed along the wire with the following distribution: 
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with ctgb the standard deviation of a GB from the uniform distribution along the wire. This distribution yields the following 
average value for the squared absolute value of the matrix element, : 
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Furthermore we have considered two extensions of this model to introduce some tilt of the GBs, either having a fixed tilt angle 
e in 


II-A or being randomly tilted in II-B 


A. GB planes tilted under fixed angle 

We consider N GBs with the same distribution as in the MS model, but with all the planes tilted under the same angle, the 
potential can be modeled by two extra parameters compared to the MS model: 
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Fig. 1. (a) The grain boundaries are all tilted w.r.t. the y-axis with the same angle 6. (b) The grain boundaries are randomly tilted, determined by the plane 

parameters jSa, 7 a, uniformly distributed in the interval [—A, +A]. 
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Fig. 2. The relaxation times (correct: blue circles, assuming equal RT: yellow squares) for all Fermi level states of a wire with cross-section of 10 by 10 
~ 3.6 nm) are shown for (a) perpendicular GB planes (6 = 90°) and (b) tilted GB planes (0 = 70°) as a function of their transport momentum 
kz, limited to positive momenta. The constant Fermi wave vector kp is shown as a red, dashed line. The simulated wire has N grain boundaries with energy 
baiTier of S = 1.5eV-acu and LzjN = D. 


The squared matrix element changes in the following way: 
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The only approximations that were made are the same as in the MS model, namely that N is quite large and ctgb ^ L^. 


30 



-■-P=0 (0-90°) 
-♦-P=0.2{e«80°) 
P=0,4 (0=70°) 
-A-P=O,6(0=6O°) 
-►-P=O.8(0=5O°) 
p=1.2(0=4O°) 


10 


45 



0123456789 

Diameter (nm) 


-♦-A=0.l 

A=0.2 

-A-A=0.3 

-►-A=0.4 

A=0.5 

10 


(a) (b) 

Fig. 3. The resistivity of simulated wires with grain boundaries planes tilted (a) with a fixed angle 9 and (b) randomly with a tilt parameter A are shown 
as a function of the wire diameter. The grain boundaries have energy barrier S = 1.5eV-acu and LzjN = D. 



















































In Fig. 1^ simulation results for resistivities of different diameter wires are shown for GB planes tilted in the cc-direction 
with several values of the tilt angle. The tilt angle 6 is in this case related to the /?o by tan6 = /3o- 


B. Randomly tilted GB planes 

If every GB plane has a different orientation, we would need two parameters per plane to describe the potential. For the 
sake of simplicity however, we consider an ensemble of random tilts Va(x,y,z) = UgbLgbS [z — — PaX — 'jau), with 

uniformly distributed Ponla G [~Ax/i//2, - The two remaining parameters that represent the GB plane tilt in x- and 

y-direction respectively are A^, and A^. Hence the tilt distribution is considered to be the same for each GB plane and uniform. 
The advantage of this distribution, as compared to a uniform distribution in angular coordinates, is that the scattering matrix 
elements can still be obtained analytically, while capturing a random tilt of each GB. With this ensemble, one has to be careful 
for large values of I^xjy however. The distribution of tilt angles will be dominated by large angles because of the Jacobian 
when transforming to spherical coordinates. 

The constants ^xjy define a maximum tilt for the GBs. The squared matrix element averaged over GB positions and 
orientations reads; 
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with Afc= V 

The resistivities obtained witn si 
Aj, are shown in Fig. (b). 
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simulations of randomly oriented GB planes for different values of tilt parameter A = A^^ = 


III. Conclusion and remarks 


We have simulated the resistivity for metal nanowires with cross sections up to lOnm by lOnm with GB planes tilted under 
a fixed angle or randomly tilted as a scattering source for the electrons. The resistivity values clearly decrease for all diameters 
when the GB planes are randomly tilted, as shown in Fig. When each GB plane is randomly tilted, there is also a lower 
resistivity w.r.t. perpendicularly oriented GB planes, but the difference is barely visible. An important aspect of these resistivity 
values is the correct retrieval of the sub-band state dependent RTs. From Fig. it is apparent that there is no clear relation 
between the transport momentum and the RT of an electron and that approximate methods to obtain the RTs can give incorrect 
results. It is very important however to have precise results for this relation, because it has a direct impact on the resistivity 
and conductivity through ([^. In the case of a fixed tilt angle for all GB planes, the approximation assuming equal RTs in Q, 
shown in Fig. (b), underestimates the conductivity, but the reverse effect, an overestimation, can also be observed if other 
scattering mechanisms are considered |j8]. In the case of tilted GB planes, a lot of very stable outliers and a stability trend 
for high transport momentum electron states are missed by the equal RT approximation, neither of these effects being present 
in case of perpendicular GB planes. Very interestingly, the RTs for high transport momentum electron states goes up and 
approaches the value obtained with another approximation method, ignoring the incoming scattering while applying Fermi’s 
golden rule. 

The two profiles of GB tilt that were presented in this text ignore some aspects of GBs. As a first remark, the resistivities 
shown in Fig. (a) only show GB planes tilted in the x-direction. When the GBs are tilted in a combination of the x- and 
y-direction, the results will be quantitatively different, but the overall trend remains the same. A second remark addresses the 
correlations between the orientations of two different GB planes. For low GB densities, it is a safe approximation to neglect 
these correlations as we have done, but for higher densities correlations between adjacent GB planes can be expected. Deviations 
from the potential form in Q can as well be modeled, definitely for smaller diameters. Further input of GB properties and 
statistics in small diameter metal nanowires is required to make more realistic GB models. Other correlations that could be 
included are those between the tilt angle and the energy barrier height and/or width. When those correlations are observed 
experimentally, it would be interesting to include them in the statistical ensemble of Vgb- If the energy barrier height or width 






















would increase substantially for tilted GBs, the decrease in resistivity for larger tilt angles could be reduced. However, it is 
expected that a more realistic description of GBs in metal nanowires would still show a resistivity decrease for increasing 
tilt angles. It is easy to show that in the limit where all GBs are tilted under 90°, we retrieve a form of surface roughness 
scattering as described by Ando’s model, which has a much lower resistivity contribution |[8). 
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